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Aostract 

The structure of fuel sprays in a combustion 
cnamoer is tneoretical ly investigated using com- 
puter models of current interest. Three repre- 
sentative spray models are considered: (1) a 

locally nomogeneous flow (LHF) model, which 
assumes infinitely fast interphase transport 
rates; (2) a deterministic separated flow (OSF) 
model, which considers finite rates of interphase 
transport out ignores effects of droplet/ 
turDulence interactions; and (3) a stochastic 
separated flow (SSF) model, which considers 
droplet/turDulence interactions using random 
sampling for turbulence properties in conjunction 
with random-walk computations for droplet motion 
and transport. Two flow conditions, i.e., swirl 
number = 0.72 and 1.98, are studied to investigate 
the influence of swirl on droplet life histories 
and the effects of droplet/turbulence interactions 
on flow properties. Comparison of computed 
results with tne experimental data show that gen- 
eral features of tne flow structure can be pre- 
dicted with reasonable accuracy using tne two 
separated-f low models. In contrast, the LHF model 
overpredicts tne rate of development of the flow. 
Wnile tne SSF model provides better agreement 
with measurements than the DSF model, definitive 
evaluation of the significance of droplet/ 
turoulence interaction is not achieved due to 
uncertainties in tne spray initial conditions. 


Nomenclature 

B mass transfer driving potential, Eq. (14) 

Cg drag coefficient 

dp droplet diameter 

D binary diffusivity 

f mixture fraction 

g mean square of mixture fraction fluctuations 

h heat transfer coefficient 

nfg heat of vaporization 

k s turbulence kinetic energy 

Le Lewis number 

mp droplet mass 

m'" mass flux at droplet surface 

n-j number flow rate of droplets in group i 

N p ,N s convection correction, Eq. (15) 

P ( f ) probability density function (POF) of f 


Pr 

Prandtl number 


r 

radial distance 


r o 

combustion chamber radius 

Re 

droplet Reynolds 

number 

S 

swirl number 


Sc 

Scnmidt number 


s* 

source term 


S Pi> 

droplet source 

term 

t 

time 


T 

temperature 


u 

axial velocity 


V 

radial velocity 
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w tangential velocity 

Vj volume of computational cell j 
x axial distance 

Yp mass fraction of fuel vapor 
r exchange coefficient 

e rate of dissipation of turbulence kinetic 

energy 

x tnermal conductivity 

u viscosity 

p density 

$ generic property 

\(i scalar property 

Supscripts 

p droplet property 

F fuel vapor 

s droplet surface property 

» amoient property 

Superscri pts 

(“) time-averaged quantity 

( ) Favre-averaged quantity 


Introduction 

The objective of this investigation is to 
evaluate liquid fuel-spray combustion models for 
gas turbine combustors and furnaces. Numerous ' 
spray models nave been proposed and some are cur- 
rently being used in industry to assist gas tur- 
bine combustor development. 1 ’ 2 Recent reviews 
oy Crowe,-* Law, Faeth 5 and Sirignano 6 
discuss important aspects of spray modeling that 
are currently being addressed. Recent spray 
models differ in details, but generally may be 
divided into two categories: tne locally homo- 

geneous flow (LHF) models and separated flow (SF) 
models. 

LHF models represent the simplest treatment 
of a multiphase flow and have been widely used to 
analyze sprays. 7 >° Tne key assumption of the 
LHF model is that interphase transport rates are 
fast in comparison to the rate of development of 
the flow as a whole. This implies that all phases 
have tne same velocity and temperature, and phase 
equilibrium is maintained at each point in the 
flow. Clearly, LHF models are only formally cor- 
rect for flows containing infinitely small drop- 
lets. The potential range of application, how- 
ever, can be much larger, since LHF models can be 
used when particle relaxation times are small in 
comparison to the characteristic time scales of 
tne continuous phase. LHF models have been 
extensively evaluated by Faeth and coworkers 8 ’ 7 
ano by Khalil and Whitelaw. 8 It was found that 
LHF models provide a reasonable description of 
tne structure of a spray, however, they generally 
overestimate the rate of development of tne flow. 

SF models have been proposed to consider 
interpnase transport phenomena. Numerous 
separated T f low models have been proposed, Gosman 
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Theoretical Methods 


and Jonns^ have divided them into three groups: 
discrete droplet models, continuous droplet models 
and continuum formulation models, l^he discrete 
droplet model was first proposed by^Crowe^ and 
later used by Gosman and Ioannides 11 and E1‘ 

Bannawy and Wnitelaw. 2 In this model droplets 
are divided into representative groups whose 
trajectories are tracked by solving the Lagrangian 
ordinary-differential equations that govern drop- 
let transport. The gas-phase conservation equa- 
tions, with source terms included to represent 
tne effect of the droplets on the gas-phase, are 
solved in tulerian coordinates. In the continuous 
droplet model, proposed by Williams 11 and used 
oy Westorook, 1 * 1 tne droplets are represented by 
a statistical distribution function in a multi- 
dimensional space of droplet size, velocity, 
location and time. The properties of droplets 
are determined by solving the conservation equa- 
tion of the distribution function, similar to the 
analysis of molecular properties in the kinetic 
theory of gases. In the continuum formulation 
models, droplets are treated as a continuous 
fluid, so that the gas and droplets are repre- 
sented as interpenetrating continua. This formu- 
lation requires one droplet fluid phase for each 
droplet size. It is very expensive in terms of 
computer storage and time for practical sprays, 
since they generally have a wide range of droplet 
sizes. In most recent spray models the discrete 
droplet formulation has been adopted, since it 
reduces numerical diffusion and computer storage 
requirements while providing a convenient frame- 
work for dealing with multiple droplet size and 
complex interpnase transport phenomena. The pre- 
sent investigation emphasizes the discrete droplet 
model . 

Many discrete droplet models neglect the 
effects of turbulence on interphase transport. 1 ’ 2 
This implies tnat droplets follow deterministic 
trajectories, yielding the deterministic separated 
flow (i)SF) model. Neglecting the effects of tur- 
Dulence on droplet transport is appropriate when 
cnaracteri stic droplet relaxation times are large 
in comparison to characteristic times of turbulent 
fluctuations. Few practical sprays, however, 
satisfy this condition. Dukowicz, 1 ^ Gosman and 
Ioannioes 11 and Faeth and coworkers 1 ®' 2 ^ have 
adopted stochastic metnods to treat the turbulence/ 
droplet interactions. In the stochastic separated 
flow (SSF) model, droplets are assumed to interact 
with turbulent eddies whose properties are- obtained 
by random sampling. An appropriate number of 
random samples is taken for each droplet group so 
that statistically significant mean droplet prop- 
erties can be obtained. The formulation of the 
SSF model was first proposed by Gosman and 
loannides 11 to study droplet dispersion. Faeth 
and coworkers 1 ^ -2 * 1 furtner developed the model 
to include tne effects of turbulence on interphase 
neat and mass transport. Tne SSF model has been 
extensively evaluated by Faeth and coworkers 1 ® -211 
in a wide variety of parabolic flows and the 
results nave been very encouraging. 

The present investigation extends the SSF 
models of Refs. 16 to 20 to consider swirling and 
recirculating flows. Results from LHF and DSF 
models are also presented in order to hignlight 
tne effects of finite interphase transport rates 
ana droplet/turbulence interactions. 


The LHF model is equivalent to the analysis 
of a variable density single-phase flow. The 
effects of the dispersed phase appear in tne 
upstream boundary conaitions and the representa- 
tion of state relationships. The present LHF 
model is very similar to that reported in earlier, 
work on evaporating ana combusting sprays . 1 9»20,22 

Gas Pnase 


The major assumptions of the gas-phase for- 
mulation are as follows: (1) the flow is steady, 

axisymmetric, and swirling, (2) exchange coeffi- 
cients of all species and thermal energy are the 
same, (3) Kinetic energy of the mean flow is neg- 
ligiDle, (4) combustion is diffusion controlled, 
(b) radiative neat transfer is neglected, and (6) 
the sprays are dilute so that droplet collisions 
and the displacement of the gas by the droplets 
can be neglected. Favre-averaged governing 
equations are solved. The conserved-scalar 
formulism is used to obtain scalar properties. 
Turbulence closure- is achieved by employing a 
k-e-g turbulence model. 12 ’ 21 

Tne governing equations of the gas-pnase can 
be written in the general form, 12 

h Gu * } + r h = h 

r 


where = ‘p^/p is the Favre-averaged dependent 
variaDle. The expressions for the exchange coef- 
ficient r and tne source term S can be found 
in Refs. 12 and 23. The term S p4 , is used for 
tne SF models and represents the interaction 
between the droplets and the gas phase ana is 
obtained from the droplet trajectory calculations, 
as will be described later. Because of the dilute 
spray assumption the influence of droplets on 
turbulence properties is neglected. 

A gaseous diffusion flame structure is 
assumed for the combustion analysis and the 
droplets are treated as distributed sources of 
fuel vapor. The presence of envelope flames 
around or wake flames behind Individual droplets 
is neglected since the majority of droplets spend 
most of their lifetimes in fuel rich regions. 
Experimental confirmation of this assumption can 
oe found in Refs. 5 and 24. 

Tne Favre-averaged scalar properties (e.g., 
temperature and species concentrations) are cal- 
culated from tne probabi 1 i tv density function 
(PDF) of mixture fraction, P(f), and state rela- 
tionships of scalar properties as a function of 
mixture fraction, ( f ) , 
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.s clipper Gaussian function is assumed for 
"S'- r j ana tne most grooaole value and variance 
are oeierminea by t ana g. State relationships 
are ootainea by thermodynamic equilibrium calcu- 
lations using tne CEC-76 computer program devel- 
oped by Gordon ana Mcbriae. 3 The SSF model 
requires tne time-averaged PDF, P(f), wnicn can 
oe found as follows: 

P(f) = p/ 0 (f) P(f) (4) 


droplet Phase 


Tne major assumptions involved in the drop- 
let transport calculations have been discussed in 
Detail Dy Law 4 and Faetn s and will not be 
repeated nere. In tne present investigation, the 
cnin-sKin approximation is usea for droplet evap- 
oration. Tnis implies that the bul< mass of a 
droplet remains at its initial condition and only 
tne surface of tne droplet is heated during evap- 
oration. Tnis approximation is convenient for 
tne conserved scalar approach of the gas field 
since energy aDSorDed Dy tne liquid is neglected. 
Tne justification for tnis approximation is dis- 
cussed in Refs. 5 ana 20. 

Droplet trajectories are solved in Lagrangian 
coordinates. The governing equations of droplet 
transport are given oelow. 

uroplet motion . The position, axial, radial 
ano tangential momentum equations are 
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N or N = 1 + 
P s 


0.276 Re 1/2 (Pr or Sc) 1/3 
Ll + 1.232/ (Re (Pr or Sc) 4/3 ]^ 2 


(15) 


Tne fuel vapor concentration and temperature at 
tne liquid surface are related through the fuel 
vapor pressure characteristics 

*Fs ■ f <V »> (16) 

Time-averaged gas properties and instantan- 
eous eddy properties should be used for droplet 
trajectory calculations in the OSF and SSF models, 
respectively. Time-averaged scalar properties 
can be found directly from the conserved-scalar 
formulation. 27 However, to convert Favre- 
averaged velocities to time-averaged velocities 
requires correlations of velocity-density fluctu- 
ations. To avoid undue complications in the 
analysis, the Favre-averaged gas velocities are 
used in the calculations. Jeng and Faeth 28 
estimated the difference between the two averaged 
velocities in their diffusion flame to be less 
than 5 percent. 

Droplet/Turbulence Interaction 


In the SSF model, the effects of turbulence 
on droplet transport are treated by a stochastic 
approach. Tne complete formulation of the stoch- 
astic method can be found in Refs. 9 and 16 to 20 
and only a brief description is given here. 


Tne standard drag coefficient is approximated as 
fo 1 lows: 
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Re < 1000 
Re > 1000 
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wnere the droplet Reynolds number, Re, is defined 


as Re = 
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P 


/ii. 


Droplet evaporation . The droplet heat and 
mass transport equations are as follows: 28 



( 10 ) 


Lagrangian equations of droplet motion and 
transport are solved by dividing the droplets 
into n groups (defined by initial position, 
size and velocity) at the nozzle exit and comput- 
ing their life histories in the flow field. The 
droplets are assumed to interact with a random 
distribution of eddies along their trajectories. 
Tne eddy properties are obtained by random sam- 
pling of the PDF’s of velocity and mixture frac- 
tion. To simplify the analysis, velocity and 
mixture fraction are assumed to be statistically 
independent. The turbulence is assumed to be 
isotropic and to have a Gaussian PDF in the 
velocity fluctuations. Although the time-averaged 
PDF should strictly be used, the Favre-averaged 
PDF is used here for the same reason as in the 
DSF analysis. 

The scalar eddy properties are obtained by 
sampling the time-averaged PDF, P(f), for an 
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instantaneous value of f and then obtaining the 
corresponding scalar properties from the state 
relationsnips . The droplet/eddy interaction time 
is assumed to be the minimum of either the eddy 
lifetime or the transit time for the droplet to 
cross the eddy. Details of the specif ication of 
tnese times and their calibration can be found in 
Refs. 16 to 20. 

Droplet Source Terms 

Tne influence of droplets on the turbulence 
is neglected in this study and no source terms 
from the droplets are included in the k, e and 
g equations. The droplet mass and momentum 
source terms appearing in Eq. (1) are computed 
after the droplet transport equations are inte- 
grated. 

Tne mass and axial momentum source terms for 
each computational cell j are calculated from 
the net change of droplet mass and axial momentum 
for each droplet group traversing tnis cell, as 
fo 1 1 ows : 


derivatives of all dependent variables are assumed 
to De zero at the exit. Axisymmetry is assumed 
and adiabatic and non-slip conditions are speci- 
fied along all solid surfaces. 

The Lagrangian formulation of the droplet 
transport equations requires the specification of 
initial conditions for the droplets to start the 
integration. The initial conditions required in 
the present analysis include aroplet velocities, 
location, and size distribution. This informa- 
tion, however, is rarely available in practical 
flows. As a result, approximations of initial 
conditions are used in the predictions. Initial 
droplet locations and velocities are estimated 
from the limited available experimental informa- 
tion. Tne droplet size is assumed to have a 
Rosin-Rainmler size distribution. Upon impingement 
on a wall, tne droplets are assumed to b ounce 
Dack elastically with an angle of reflection equal 
to tne angle of incidence. 


Numerical Solution Procedure 



The droplet radial and tangential momentum source 
terms, however, cannot be calculated in this man- 
ner due td the additional force terms in the 
respective droplet equations of motion (cf. Eqs. 
(7) and (8)). These source terms are obtained by 
directly computing the change of radial and tan- 
gential momentum due to drag, as follows: 



A modified version of the TEACH computer 
program is used to solve the gas-phase governing 
equations. It includes an improved finite- 
difference procedure, the Bounded Skew-Upwind 
Differencing (BSUD) method developed Dy Raithby.30 
Tne pressure field is estimated using the Pressure 
Implicit Split Operator (PI SO) predictor-corrector 
technique developed by Issa. The axisymmetric 
flow domain is divided into relatively fine 
meshes, i.e., 60 by 60. For the separated-f low 
analysis, the ordinary differential equations 
governing the motion and transport of each droplet 
group are integrated using a fourth-order Runge- 
Kutta method. The numerical integration of drop- 
let trajectories is started from the fuel nozzle 
exit where the initial conditions of the droplets 
are prescribed. Seven hundred and 2800 trajec- 
tories are computed for the DSF and SSF models, 
respectively, in all cases. The solution proce- 
dure follows the PSI-CELL (Particle-source in 
Cell) iterative metnoo of Crowe et al. 10 The 
gas-pnase solution is oDtained with droplet source 
terms from the previous iteration, then the drop- 
let trajectories and tne new mass and momentum 
source terms in each cell are calculated. The 
iteration cycle is repeated until the solutions 
for both phases converge. Generally, convergence 
is achieved in 1200 interations for the separated- 
flow calculation and 750 iterations for the LHF 
calculation. The CPU time requirements are 
approximately 600, 1400, and 2800 sec, respec- 
tively, for tne LHF, DSF, and SSF models on the 
CRAY-1S/2200 computer. 


X (w . - w) 

u . - u 

1 ) 

pi 

pi 

1 / 


Boundary and Initial Condition s 

Tne elliptic nature of the gas-phase govern- 
ing equations requires that boundary conditions 
be specified along all surfaces enclosing the 
computational domain. Measured values of mean 
axial velocity and temperature are prescribed at 
tne inlet. Inlet radial velocity is taken as 
zero and the tangential velocity is specified to 
have a solid-body rotation profile. Second 


Results and Discussion 

Tne theoretical models described previously 
were applied to the liquid-fueled furnace used by 
Knalil et al.“ The geometry of the combustion 
cnarrtDer is shown in Fig. 1. Two swirl numbers 
(S = 0.72 and 1.98) were considered and the fuel- 
air equivalence ratio was fixed at 0.8. 

The predicted and measured temperature pro- 
files at different axial locations are given in 
Fig. 2 for S = 0.72. The results snow that both 
the OSF and SSF models underpredict the temperature 
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in the near inlet region, while the LHF model 
overpredicts the temperature, especially at large 
radii (r/r 0 > 0.8). The discrepancies between 
tne predictions of the two separated-f low models 
and the measurements are due to the underestima- 
tion of droplet spreading and evaporation, which 
is probably caused by uncertainties in the spray 
initial conditions. The SSF model snows better 
agreement than the DSF model because of better 
predictions of droplet evaporation and spreading 
rates. The overestimation of temperature by the 
LHF model is due to tne assumption of infinitely 
fast interphase transport. 

Far downstream of the combustor inlet, all 
three models predicted higher temperatures than 
tne measurements. This is attributed to the 
neglect of radiative heat transfer. However, 
these errors in the predicted temperatures do not 
significantly affect the performance of the three 
models and hence, the conclusions of this study. 
Both tne experimental data and the separated-f low 
model predictions showed high temperature peaks 
near the combustion chamber wall at downstream 
locations. The peak temperatures occur in regions 
where droplets have traversed and evaporated. The 
LHF model did not predict a temperature peak near 
the chamber wall because this model does not con- 
sider discrete droplet motion. 

Figure 3 presents the experimental data and 
predictions for S = 1.98. Tne temperature pro- 
files are more uniform for S = 1.98 than for 
S = 0.72 (Figs. 2 and 3). This is due to the 
increased mixing obtained with increased swirl. 

The nigner swirl also produces a stronger recir- 
culating flow (i.e., more fuel vapor is transpor- 
ted upstream) and hence higher temperatures are 
obtained in the upstream region (x < 0.24 m) for 
S = 1.98. The agreement between the predictions 
and tne measurements is similar to the lower swirl 
number case, except close to the inlet at large 
radii, where both the DSF and SSF models signifi- 
cantly underestimate the temperature. Very few 
fuel droplet groups are predicted to arrive in 
these regions (i.e., close to the inlet at large 
raoii). Tnis is probably caused by inaccurate 
initial droplet conditions. The fuel vapor con- 
centration in these regions is therefore under- 
predicted which results in much lower temperatures 
than the experimental data. For both swirl number 
cases in regions far downstream, the DSF model 
predicts much higher temperature peaks, less uni- 
form temperature distributions, and larger dis- 
crepancies with measurements than the SSF model 
(Figs. 2 and 3). This is because the DSF model 
neglects turbulent dispersion and hence, under- 
estimates droplet spreading. 

Radial profiles of gas-phase axial velocity 
are presented in Figs. 4 and 5 for S = 0.72 and 
1.98, respectively. The boundaries of the central 
recirculation zones for S = 1.98, 1.25 and 0.94 
are shown in Fig. 6. For clarity in presentation, 
the results of only the SSF model are given. It 
can be seen that Doth tne magnitude of the reverse 
velocities and the size of the recirculation zones 
are underpredictea, the errors increase with 
increasing swirl number. The implication of the 
underestimation of the recirculation can be seen 
in the temperature distributions in Fig. 3. Tne 
experimental data shows more uniform distributions 
oecause of greater mixing caused by larger reverse 


velocities. Tne poor performance of the model in 
predicting the size of the central recirculation 
zone and the reverse velocity may be partly attri- 
buted to the k - e model used in the present 
analysis. It is well known that k - e models 
nave not been very successful in predicting highly 
swirling flows. ” Beyond the central recircu- 
lation zone, the predicted axial velocities are 
higher than the measured values due to the higher 
predicteo temperatures. 

Tne predicted and measured radial profiles 
of gas-phase tangential velocities for the higher 
swirl numDer case are illustrated in Fig. 7. The 
calculated tangential velocities are smaller than 
tne measured values at almost all locations. 

This discrepancy is attributed to uncertainties 
in the inlet swirl velocities. It should also be 
noted that the swirl numbers quoted by Knalil 
et al.^ for the experimental data were esti- 
mated from the swirler vane angle rather than 
direct measurements and may introduce some error. 
The underestimation of tne tangential velocities 
is partially responsible for the errors in the 
predictions of flow reversal. 

The predicted axial profiles of liquid mass 
flow rate are shown in Fig. 8 for S = 0.72 and 
1.98. As expected, the LHF model produces unre- 
alistically high evaporation rates. This 
high rate of evaporation is due to the method 
used to obtain the inlet boundary conditions and 
the assumption of infinitely fast interphase tran- 
sport rates. The results of the two separated- 
flow models appear reasonable, but lack of 
experimental data prevents quantitative compari- 
sons. The SSF model predicts higher evaporation 
rates than tne DSF model Decause it considers the 
effects of turbulence on interpnase transport, 
wnich enhances the rate of evaporation. 20 
Tne influence of swirl strength on the overall 
evaporation rate is clearly indicated in Fig. 8, 
which shows faster evaporation for higher inlet 
swirl. Tnis can oe explained by the interaction 
Detween the central recirculation zone and the 
droplet trajectories. In the high swirl- case a 
large reverse flow is generated and most small 
droplets are trapped in the recirculation zone 
and evaporate there, resulting in intensive burn- 
ing and high temperatures in tnis region. Large 
droplets which penetrate the central recirculation 
zone burn further downstream or near the combustor 
wall. In the low swirl case, most droplets sur- 
vive the central recirculation zone and continue 
to evaporate in the far downstream region. Since 
the mixing is weak outside the reverse flow 
region, high fuel vapor concentration along the 
droplet trajectories produce local high tempera- 
ture peaks (see Fig. 2). In addition to the 
influence on droplet trajectories, the larger 
reverse flow zone recirculates more hot combustion 
products from downstream regions and increases 
tne temperature at near inlet locations. This in 
turn contributes to a higher evaporation rate for 
the hign swirl case. 

The predicted axial profiles of droplet SMD 
are illustrated in Fig. 9 for S = 0.72. The 
values presented here were obtained by integrating 
over the cross-sectional area at each axial loca- 
tion. Following Khalil et al.,“ the spray was 
assumed to have an initial SMD of 127 urn. A 
Rosin-Rammler size distribution, with 15 size 
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groups and size range of 10 to 290 ym, was 
employed in the calculations. As the droplets 
move downstream, the SMO increases and then 
aecreases. The initial increase of SMD is caused 
Dy the reduction of the number of small droplets, 
due to the preferential evaporation of small drop- 
lets (i.e., small droplets evaporate faster than 
large droplets). Wnen most of the small droplets 
have completely evaporated, the SMD begins to 
decrease because of the continuous evaporation of 
large droplets. The DSF and SSF models produce 
similar results, but tne SSF model produces more 
rapid variations in the SMD than the DSF model 
Decause it predicts a higher evaporation rate. 


Conclusions 

Current theoretical spray models have been 
used to predict properties of fuel sprays in a 
cylindrical combustion chamber. The results 
showed that general features of the gas-phase 
flow field can be predicted reasonably well. 
However, the evaluation of model performance for 
the dispersed phase is not conclusive mainly due 
to inadequate knowledge of the droplet initial 
conditions and a lack of experimental data for 
downstream droplet properties. Some quantitative 
differences between predictions and me 'surements, 
such as size of the recirculation zone, magnitude 
of the tangential velocity, temperature distribu- 
tion close to the inlet, could be related to 
uncertainties in the gas-phase upstream boundary 
conditions and also to droplet initial conditions. 
Tne need for a set of benchmark experimental data 
on comDustor fuel sprays with characterization of 
inlet conditions and flow measurements in both 
pnases has long been recognized. 3,11,12 j^e 
present investigation reemphasizes this need. 

The LHF method yields unrealistically high 
evaporation rates and overestimates the develop- 
ment of the flow field. Therefore, it has only 
limited utility in modeling practical sprays. 

The dispersion model gives less concentrated drop- 
let trajectories and consequentially better temp- 
erature distributions than tne deterministic 
metnod. Turbulence/droplet interaction enhances 
the rate of evaporation and produces higher fuel 
vapor concentration in the upstream portion of 
the combustor, where most droplets evaporate. 
Inclusion of these effects results in better 
predictions of the temperature, especially in 
regions near the wall and close to the inlet. 
Although the SSF model appears to have certain 
advantages over DSF model, the differences 
between the two predictions are not large 
compared to tne uncertainties in the droplet 
initial conditions. The lack of experimental 
data providing important inlet conditions and 
detailed flow measurements in both phases 
prevents a decisive evaluation of the two 
separated-f low models. 


Acknowledgement 

This research was sponsored by the National 
Aeronautics and Space Administration, Contract 
No. NAS 3-24105, under the technical management 
of Daniel L. Bulzan of the Lewis Research Center. 


1. Swithenbank, J., Turan, A., and Felton, P.G., 
Das Turbine Combustor Design Problems, 
Lefebvre, A.H. ed., pd. 249-314. Hemisphere 
Publishing, Washington, 1980. 

2. Mongia, H.C. and Smith, K. “An Empirical 
Analytical Design Methodology for Gas Turbine 
Combustors." A1AA Paper No. 78-988, 1978. 

3. Crowe, C.T., “Review-Numerical Models for 
Dilute Gas-Particle Flows," J. Fluid Enqrq., 
Vol . 104, 1982, pp. 297-503. 

4. Law, C.K., "Recent Advances in Droplet 
Vaporization and Combustion," Prog. Energy 
Combust. Sc. , Vol. 8, 1982, pp. 171-201. 

5. Faeth, G.M., "Evaporation and Combustion of 
Sprays," Prog. Energy Combust. Sci., Vol. 9, 
1983, pp. 1-76. 

6. Sirignano, W.A., "Fuel Vaporization and Spray 
Combustion Theory," Prog. Energy Combust. 
Sci. , Vol. 9, 1983, pp. 291-322. 

7. Shearer, A.J., Tamura, H., and Faeth, G.M., 
"Evaluation of a Locally Homogeneous Flow 
Model of Spray Evaporation," J. of Energy, 
Vol. 3, 1979, pp. 271-278. 

8. Khalil, E.E. and Whitelaw, J.H., "Aerodynamic 
and Thermodynamic Characteristics of 
Kerosene-Spray Flames," Sixteenth Symposium 
(International) on Combustion , Tne Combustion 
Institute, Pittsburgn, 1977, pp. 569-576. 

9. Gosman, A.D. and Johns, R.J.R., "Computer 
Analysis of Fuel-Air Mixing in Direct- 
Injection Engines," SAE Paper No. 800091, 
1980. 

10. Crowe, C.T., and Stock, D.E., "A Computer 
Solution for Two-Dimensional Fluid-Particle 
Flows," Int. J. Num. Method. Engr., Vol. 10, 
1976, pp. 185-196. 

11. Gosman, A.D. and‘ Ioannides, E r, "Aspects of 
Computer Simulation of Liquid-Fueled 
Combustors," AIAA Paper No. 81-0323, 1981. 

12. El Banhawy, Y., and Whitelaw, J.M., 
"Calculation of the Flow Properties of a 
Confined Kerosene-Spray Flame, - " AIAA J. , 

Vol. 18, 1980, pp. 1503-1510. 

13. Williams, F.A., Combustion Theory , Addision- 
Wesley, Reading, MA, 1965. 

14. Westorook, C.K., "Three Dimensional Numerical 
Modeling of Liquid Fuel Sprays," Sixteenth 
Symposium (International) on Combustion , The 
Combustion Institute, Pittsburgh, 1977, 

pp. 1517-1526. 

15. Dukowicz, J.K., "A Particle-Fluid Numerical 
Model for Liquid Sprays," J. Comp. Phys., 

Vol. 35, 1980, pp. 229-253. 


6 



16. Shuen, J.S., Solomon, A.S.P, Zhang, Q.F. and 

Faeth, G.M., "Structure of Particle-La'den 
Jets: Predictions and Measurements," AIAA J ■ , 

Vol . 23, 1985, pp. 396-404. 

17. Solomon, A.S.P. , Shuen, J.S., Zhang, Q.F., 

and Faetn, G.M., "Structure of Nonevaporating 
Sprays, Part I.: Initial Conditions and Mean 

Properties," AIAA J. , Vol. 23, 1985, 

pp. 1548-1555. 

18. loio., "Structure of Nonevaporating Sprays, 

Part II: Drop and Turbulence Properties," 

AIAA J. , Vol. 23, 1985, pp. 1724-1730. 

19. iDio., "Measurements and Predictions of the 
Structure of Evaporating Sprays, J. Heat 
Transfer , Vol. 107, 1985, pp. 679-686. 

20. Shuen, J.S., Solomon, A.S.P., and Faeth, G.M., 
"Drop-Turbulence Interactions in a Diffusion 
Flame," AIAA J. , in press. 

21. Lockwood, F.C. and Naguib, A.S., "The 
Prediction of the Fluctuations in the 
Properties of Free, Round Jet, Turbulent, 
Diffusion Flames," Comb. Flame, Vol. 24, 

1975, pp. 109-124. 

22. Mao, C.P., Szekely, G.A., Jr., and 
Faeth, G.M., "Evaluation of a Locally 
Homogeneous Flow Model of Spray Combustion," 

J. Energy , Vol. 3, 1980, pp. 79-87. 

23. Gosman, A.D., loannides, E., Lever, D.A., and 
Cliffe, K.A., "A Comparison of Continuum and 
Discrete Droplet Finite-Difference Models 
Used in tne Calculation of Spray Combustion 
in Swirling Turbulent Flows," AERE Harwell 
Report TP 865, 1980. 

24. Onuma, Y. and Ogasawara, M., "Further 
Experiments on a Spray Combustion Flame," 
Sixteenth Symposium (International) on 
Combustion , The Combustion Institute, 
Pittsburgn, 1977, pp. 561-576. 


25. Gordon, S. and McBride, B.J., "Computer 
Program for Calculation of Complex Chemical 
Equilibrium Compositions, Rocket Performance, 
Incident and Reflected Shock, and 
Cnapman-Jouget Detonations," NASA SP-273, 

1976. 

26. Faeth, G.M., "Current Status of Droplet and 
Liquid Combustion," Proq. Enerqy Combust. 

Sci. , Vol. 3, 1977, pp. 191-224. 

27. Bilger, R.W., "Turbulent Jet Diffusion 
Flames," Proq. Enerqy Combust Sci., Vol. 1, 
1976, pp. 87-109. 

28. Jeng, S.M. and Faeth, G.M., "Species 
Concentrations and Turbulence Properties in 
Buoyant Methane Diffusion Flames," J. Heat 
Transfer , Vol. 106, 1984, pp. 721-727"! 

29. Knalil, K.H., El Mahallawy, and Moneib, H.A., 
"Effect of Combustion Air Swirl on the Flow 
Pattern in a Cylindrical Oil Fired Furnace," 
Sixteenth Symposium (International) on 
Combustion , The Combustion Institute, 

Pi ttsDurgn, 1977, pp. 135-143. 

30. Raithby, G.D., "Skew-Upstream Differencing 
Scneme for Problems Involving Fluid Flow," 
Computer Methods in Applied Mechanics and 
Engineering' ^ Vo 1 . 9, No. 2, 1976, pp. 153-164 . 

31. Issa, R.I., "Numerical Methods for Two- and 
Three-Dimensional Recirculation Flows," 
Computational Methods for Turbulent, 

Transonic, and Viscous Flows, Essers, J.A. 
ed.. Hemisphere Publishing, New York, 1983, 
pp. 183-212. 

32. Hutchinson, P., Khalil, E.E., and 
Whitelaw, J.H., "The Measurements and 
Calculation of Furnace-Flow Properties," 

J. Enerqy , Vol. 1, 1977, p. 212. 


7 






Figure 3.- Radial profiles of temperature, 
S = 1.98. 



.8 


.6 


.4 


.2 


r/r 


o 


lure 4. - Radiaf profiles of axial velocity, S ■ 










210 r~ 



0 .1 .2 .3 .4 .5 .6 


X,m 

Figure 9. - Axial profiles of Sauter mean diameter. 




1. Report No. NASA CR-1 75028 . 
AIAA-86-0450 

2. Government Accession No. 

3. Recipient’s Catalog No. 

4. Title and Subtitle 

5. Report Date 

Prediction of the Structure 

of Fuel Spravs 

December 1985 

In Gas Turbine Combustors 

6. Performing Organization Code 

7. Authors) 

8. Performing Organization Report No. 

Jlan-Shun Shuen 

None 


10. Work Unit No. 

9. Performing Organization Name and Address 

’ 


11. Contract or Grant No. 

Sverdrup Technology, Inc. 


Lewis Research Center 

NAS 3-24105 

Cleveland, Ohio 44135 

13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Contractor Report 

National Aeronautics and Space Administration 

14. Sponsoring Agency Code 

Washington, D.C. 20546 



505-31-04 ( E — 2 8 3 7 ) 


15. Supplementary Notes 


Final report. Project Manager, Daniel L. Bulzan, Internal Fluid Mechanics 
Division, NASA Lewis Research Center, Cleveland, Ohio 44135. Prepared for the 
24th Aerospace Sciences Meeting, sponsored by the American Institute of 
Aeronautics and Astronautics, Reno, Nevada, January 6-9, 1986. 

16. Abstract 

The structure of fuel sprays In a combustion chamber Is theoretically Investi- 
gated using computer models of current Interest. Three representative spray 
models are considered: (1) a locally homogeneous flow ( LHF) model, which assumes 

Infinitely fast Interphase transport rates; (2) a deterministic separated flow 
(DSF) model, which considers finite rates of Interphase transport but Ignores 
effects of droplet/turbulence Interactions; and (3) a stochastic separated flow 
(SSF) model, which considers droplet/turbulence Interactions using random sam- 
pling for turbulence properties In conjunction with random-walk computations for 
droplet motion and transport. Two flow conditions are studied to Investigate 
the Influence of swirl on droplet life histories and the effects of droplet/ 
turbulence Interactions on flow properties. Comparison of computed results with 
the experimental data show that general features of the flow structure can be 
predicted with reasonable accuracy using the two separated flow models. In con- 
trast, the LHF model overpredicts the rate of development of the flow. While the 
SSF model provides better agreement with measurements than the DSF model, defini- 
tive evaluation of the^slgnlf Icahce of droplet/turbulence Interaction Is not 
achieved due to uncertainties In the spray Initial conditions. 


17. Key Words (Suggested by Authors)) 18. Distribution Statement 

Fuel sprays; Swirling and recirculating Unclassified - unlimited 

flows; k-e turbulence model; 2-phase STAR Category 07 

f 1 ows 

19. Security Classif. (of this report) 20. Security Classif. (of this page) 21. No. of pages 22. Price* 

Unclassified Unclassified 

‘For sale by the National Technical Information Service, Springfield, Virginia 22161 



















National Aeronautics and second class mail 

Space Administration 

Lewis Research Center address correction requested 

Cleveland. Ohio 44135 



Official Business 
Penalty for Private Use $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 


NASA 




